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(<~*) ■ We investigate quantum dynamics in the electron-nuclei coupled spin system in quantum dots 

' and clarify the fundamental features of quantum correlation induced via successive electron spin 

, measurements. This quantum correlation leads to interesting phenomena such as the bunching 

of outcomes in the electron spin measurements and the revival of an arbitrary initial electron spin 
state. The nuclear spin system is also affected by the quantum correlation and is in fact squeezed via 
conditional measurements or postselection. This squeezing is confirmed by calculating the increase 
in the purity of the nuclear spin system. Thus the successive electron spin measurements provide 
a probabilistic method to squeeze the nuclear spin system. These new features are predicted not 
only for the case of a double quantum dots occupied by a pair of electrons but also for the case of 
a single quantum dot occupied by a single electron or a pair of electrons. 

■ PACS numbers: 73.21. La, 71.70.Jp, 76.70.-r, 03.67.Pp 
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I. INTRODUCTION 



Quantum state control in solid systems is a challenging task due to the strong coupling of solid state systems to 
d ' environments in contrast with the atomic systems in which the coupling to environments is much weaker. However, the 
prospect of realizing scalable architectures for quantum information processing motivated the investigation on solid 
state/semiconductor structures. Electron spins in semiconductor quantum dots(QD) proved to be one of the most 
' O ■ promising two-level systems for the quantum state control^ due to their long decoherence times. Main decoherence 
^ , mechanisms are the coupling to phonons via the spin-orbit interaction and the hyperfine(HF) interaction with the 
host nuclei. The spin-orbit interaction leads to an exponential decay of the longitudinal and transverse electron spin 
components characterized by Ti and T2 times^'^ . Under the strong confinement and a weak magnetic field, the phonon- 
mediated decoherence is greatly suppressed, whose time constant was demonstrated to reach up to lOOmsii. Instead 
. ^ the contact hyperfine(HF) interaction of the electron spin with the lattice nuclei dominates the decoherence^i^ii. 
. Contrary to the spin-orbit-mediated decoherence, the HF interaction can lead to the pure dephasing and it features 
' a Gaussian decay. The HF interaction acts on a time scale proportional to the square root of the number of nuclei 
J — , fi/T2 = A/^/N, where A is the material specific HF coupling constant, and N is the number of host nuclei. For 
T— I • example, for GaAs A = 90/.tcV^ and for a QD having 10^ nuclei, the HF induced decoherence time T2 is ~ 10ns. In 
' , order to suppress the HF induced decoherence, there have been made several proposals, such as the measurement of 
■ the HF field^ d°d^ and the polarization of nuclear spins, which will reduce the fluctuations in the HF fieldi^. However, 
00 \ to achieve these, one has to do highly precise measurements or to polarize the nuclear spins to a high degree. 

. In small mesoscopic structures the HF interaction is so far the only mechanism to probe nuclear spins, since typically 
' NMR signals from such small ensembles of nuclear spins is too weak to detect. Coherent manipulation of mesoscopic 
• ^ , ensemble of nuclear spins has been realized in semiconductor point contact devices where magnetization of nuclear 
' spins is probed by resistance measurements^^. Hyperfine interactions lead to many interesting effects, such as lifting 
^ of spin blockade in transport through double QDs^^, oscillatory currents in the spin blockade regime driven by the 
■ - - ■ HF field^^, and probing nuclear spin relaxation in the Coulomb blockade regime^^. Coherent manipulation of the spin 
state of a pair of electrons on a double QD has been achieved via electrical control of the exchange energy difference 
between the singlet and triplet spin states^, where the singlet-triplet mixing via the HF interaction has been observed. 

In light of these recent progress in the studies on the electron spin qubits and the HF interaction in QDs, we are 
going to investigate the quantum dynamics of the electron-nuclei coupled spin system, especially the manipulation and 
preparation of nuclear spin states via the HF interaction, which in turn lead to interesting effects, such as bunching 
in electron spin measurements and the electron state reviva l^^d^ . 

Our paper is organized as follows. In Sec. II, we are going to discuss a double QD model and the HF interaction 
and make comparison with available experimental data to derive relevant physical parameters. In Sec. Ill, we study 
the bunching in electron spin measurements, which arises as a result of correlations between the successive electron 
spin measurements induced by the HF interaction. In Sec. IV, we will show that nuclear spins can be conditionally 
purified via electron spin measurements, the manifestation of which is the revival of the electron spin state, enabling 
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the retrieval of an arbitrary electron spin state. These newly predicted phenomena, bunching and revival, are not 
necessarily restricted to the case of an electron pair in a double QD and can be observed in more general cases. 
In Sec. V, we discuss the feasibility to observe these phenomena in the electron spin measurements for a single QD 
occupied by either a single electron or a pair of electrons. Finally, our results and predictions are summarized in 
Sec. VI. 



II. DOUBLE QUANTUM DOT MODEL 



We are going to consider a laterally coupled double QD system occupied by two electrons. QDs are formed on a 
two-dimensional electron gas under a uniform magnetic field and the dynamics is assumed to take place only in the 
transverse spatial coordinates denoted by x and y. In this section Zeeman energies are not taken into account because 
they are not essential for the orbital dynamics. The orbital motion of electrons are governed by the Hamiltonia n^°i^^ : 

Y^{iPll^t^ + V{x,^V^,)} + V^{\p,-p,\) (1) 

i=l,2 

V{X,V) = ^^(^' - + \^^'y' - £^ - ydr) = ^ , (2) 

A(r) = |(-y,.T,0) (3) 

where the confining potential is modeled by a double well potential which can be approximated by a harmonic potential 
near x = ±a, the two-dimensional vector is represented by p= {x,y), e the external electric field , e the elementary 
electric charge (e > 0), c the light velocity in vacuum and k is the dielectric constant. Assuming the low temperature 
such that Huj kT, we study the dynamics within the manifold of the ground state orbitals, consisting of |20, S"), 
|11, To,±), |11, S) and |02, S). Here |nm, S{T±fi)) denotes the state with the electron occupation number n{m) in the 
left(right) dot and S and T±^o indicate respectively the singlet and triplet spin states. When both electrons are in the 
same dot, they are always in the singlet state since the orbital part is symmetric. However, when they are in different 
dots, the orbital part may be in an antisymmetric or a symmetric combination of ground state orbitals of the left 
and right QDs and thus the spin state may be a triplet or singlet state. For example, the orbital part of S{T)) is 
given by 

[M^)M2)±M2)M'^)]/V2 (4) 

for the symmetric(-|-)/antisymmetric(— ) combination of two electrons: one electron is localized in the left and the 
other in the right QD. The orbital state (2) is an eigenstate of the Hamiltonian |T]) approximated by the local 

harmonic potentials excluding the Coulomb potential: 



:A(rO) 
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V — — ^ '-^ + Vo{xi+a,yi) + Vo{x2-a,y2)-e{xi+X2) (5) 

4=1,2 

with Vo{x,y) = ima;^(x^ +2/^)- (6) 
Then the ground eigenstate is given by displaced harmonic oscillator states : 

0L(l)</'i?,(2) = + a - ^)0(yi)0(a:2 -a- -^)(t>{y2) , (7) 



where each wavefunction satisfies 



^^^1"^^"^^ + ima;2((x ± + y')) - ex\ ± a - ^ 
Am I muj'^ 



^(^nn±ea^-^)q^ix±a^^)<Piy) , (8) 
1 {x - xoY 

with (j}{x - xq) = — p= exp[ — ] exp[-iea;oSy/ {2hc)] , (9) 
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For the orbital parts of |02. S) and \20, S) one has to calculate the eigenstates of two electrons occupying a single 
QD including the Coulomb potential: 



«i Pi - P21 



(11) 



When the onsite Coulomb energy is smaller than the orbital energy splitting, the |02, S*) or |20, S") state orbital can 
be approximated by the product of the ground state orbitals of the harmonic oscillator. In this case, the ground state 
energy for the \02,S) and |20, S") state is given by 



2 2 
Ea2 = 2hn - 2£a ^ + 5c , E20 = 2MI + 2ea —„ 

where the onsite Coulomb energy 5c is calculated as 

5c = ^f!r/2e^/{Kl). 



+ 5c , 



The exchange energy is found by calculating the energy of the |11, S{T)) state using the full Hamiltonian ([T|): 



with 



(11, S{T)\H\ll, S{T)) ^hn + Ec±Ex± (toRiM^L) + toUMM) 

Ec = (0fl,(l)'/'L(2)|^ 1 10,^(1)^^(2)) , 

K\ri - r2\ 

(0fl(l)0L(2)|^^ AM2)M^)) , 



E 



X — 

toR — 



K\ri - r2\ 



1 Tfldy^ 



2 ^2^2 



1 



Y mw^(a; + ay , 



(12) 

(13) 

(14) 
(15) 

(16) 

(17) 

(18) 



where +(— ) corresponds to the |S')(|r)) state, Ec is the direct Coulomb energy, Ex the exchange integral and ton is 
almost equal to toL when a 3> ejimuP') and to defined by to = ^oi? — ^ol has a meaning of the single particle tunneling 
amplitude. Restricting the Hamiltonian to the relevant two electron states, it is given as 

n =2hm + {-2ea + Sc)\02, S) (02, S\ + {2ea + Sc)\20, S) (20, S\ 



-{Ec~b E |ll,T^.>(ll,r,| + (i?c + ^)|ll,^)(ll,5| 



cr=±,0 



with 



+ti?X|ll, S')(02, S\ + h.c.) + tidll, S')(20, S\ + h.c.) 
tR^{ll,S\5VR{xi) + 5VR{x2)\m,S) , j^2Ex+Ua{cpR\cpL) 
tL = {n,S\dVL{xi) + dVLix2)\20,S), 



(19) 



(20) 



where tR and are the tunneling amplitudes and it can be shown that Ir :^ ti, ^ t = \/2to, when the onsite Coulomb 
energy is smaller than the orbital energy splitting and the |02, 5) or |20, S) state orbital reduces to the product of the 
ground state orbitals of the harmonic oscillator. The structure of this Hamiltonian can be seen clearly in the matrix 
form: 



/ |20,5) |02,5) 111, 5) |ll,r,) \ 

2ea + 5c t ' 

-2ea + 5c t 

t t Ec+j/2 

V Ec -j/2/ 



(21) 



When the energy offset between the two QDs by the electrical bias is quite large, namely, ea 3> \t\, one can consider 
the dynamics only in the (1, 1) and (0, 2) charge states, where the energy of the state |11, S) is renormalized by 



5E = t^/{-2ea + Ec + j/2 - 5c) 
in consequence of the adiabatic elimination of the (2, 0) charge state. 



(22) 
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The Hamiltonian (jl9p can be put in a simpler form: 

H = -A/2|02, S){Q2, S\ + A/2|ll, S){n, S\ + t(\n, S){Q2, S\ + h.c.) 

+ {A/2- J -5E) V |ll,T,)(ll,r,|, (23) 



cr=0,± 



with A^2ea + Ec-Sc+ j/2 + 6E , (24) 

which is offset by some constant energy with respect to (|T9|) . |02, S") and |11, S*) charge states hybridize to form new 
eigenstates \ — ,S) and \ + ,S) given by 

\±,S) = ^ - (A/2T VAV4 + t2)|02,^)] (25) 



and the Hamihonian is rewritten as 



H=J^+f^[\+,S){+,S\-\-,S){-,S\]+i^-j-SE) J2 |ll,r.)(ll,r,|. (26) 

When A > \t\, \ + i~),S) |11(02),S'), whereas when A is negative and |A| > \t\, \ + (-),£') ^ |02(11),S'). The 
energy difference between the singlet ground state and the triplet states is 



J = A/2 - j - 5^; + v/AV4 + t2 (27) 

and this energy will be called the "exchange energy" in the following. In this expression "j + 6E" term coming from 
the bare exchange integral and the level shift due to the transfer integral appears and its magnitude will be estimated 
from the comparison of J with experimental data. For vanishing external magnetic field the exchange energy should 
always be positive^, namely the ground state is always a singlet state. However, in the presence of a magnetic field, 
a singlet-triplet crossing takes place at some particular value of the magnetic field, yielding a triplet ground state^^, 
i.e., J < 0. 

A. Hyperfine Interaction 

Now we are going to discuss the effects of the HF interaction with nuclei. The HF interaction is mainly described 
by the Fermi contact interaction^^: 

Ffff = A«o^S, ^.^(r.-R^). (28) 

Here ri denotes the position of the ith electron, and Rq is the position of the nucleus a. A is a material specific 
coupling constant and for instance for GaAs A = 90/ieV and vq is the unit cell volume. S and / are the spin angular 
momenta of the electron and the nucleus, respectively. When two electrons are in the same QD, they experience the 
same HF field, which implies a vanishing HF field for singlet states(which is not the case for triplet states). On the 
other hand, when the electrons are in different QDs, the mean HF field induces mixing within triplet states, and the 
difference of the HF fields in two QDs induces coupling between the singlet and triplet states. For two electrons in 
|ll,7o.±) and IlliS") states, the HF interaction is given as 

Vhf =\{hL + hn) ■ (Si + S2) + \{hL - hji) ■ (Si - S2) (29) 

= :^(|ll,T+)(ll,T+|-|ll,r_)(ll,T_|) 

^ =(/i_|ii,T+)(ii,To| + ;i_|ii,ro)(ii,r_| + h.c.) 



2V2 

^ =(-5/i_|ll,T+)(ll,5| +(5/H|ll,r„)(ll,5|+h.c.) 



2^/2 

with h = hi + hfl , 6h = hL hji, h^j^) = A vo^\(l)L{R){Ra)\'^'ta , (31) 



^idK\n,S){U,To\+h.c.) (30) 
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where ii-L(R) is the HF field in the left(right) QD and has the dimension of energy. Thus the HF fields h and Sh 
also have the dimension of energy. In general the nuclear Zeeman energy is very small, for example, for ^^Ga with 
gN — 2.02 it is about 0.74mK at i? = IT. Thus, for higher temperatures nuclear spins are randomly oriented and the 
HF field features a Gaussian distribution with the mean square value: 

K^R)) = ^0 E \'I^LiR){Ro.ri{I + 1), (32) 

where / is the magnitude of the nuclear spin and (...) denotes the ensemble average. In particular, for the uniform 
coupling, we have 

,/{hl^) - A V/(/ + l)/.yA^ , (33) 

where iV^(jj) is the number of nuclear spins in the left(right) dot. 

When the electron Zeeman energy is much larger than the HF fields, the coupling terms among the triplet states 
and those between T± and the singlet state S can be neglected and the HF interaction reduces to 

Vhf = 6h,{\n, S){n,To\ + h.c.)/2, (34) 

with Shz = hLz — hiiz being the difference of the HF fields along the applied field direction. All other spin states are 
unaffected by the HF interaction. 

Two electron system on a double QD is initialized in |02, S) state under the condition that A ^ |t|. If the bias 
voltage is changed adiabatically so that the singlet state remains in the ground state | — ,5') without ever populating 
|+, the double QD electronic Hamiltonian including the HF interaction is cast into the form: 

H + Vhf = ^(|-,5)(ll,To| +h.c.) + ^(|ll,ro)(ll,ro| - |-,5)(-,5|) (35) 
= JSz + rShzSx (36) 
with r = t/^J (A/2 + + t2)2 + ^2 ^ (37) 

where the factor r = (11, S\~, S) determines the HF coupling strength of the singlet ground state (|25|) to the triplet 
state. Now we examine the limiting values of r. The parameter A can be controlled by the bias voltage through e in 
([M)) and t can be varied through the spatial overlap of wavefunctions. When A is positive and A 3> \t\, two electrons 
are almost localized in the right dot forming a spin singlet pair and r — > 0. On the other hand, when A is negative 
and \t\ ^ |A|, two electrons are separated in different dots with negligible spatial overlap. Then \r\ — > 1 and the HF 
interaction is maximized. In Eq. ([36]) the Hamiltonian is written in the pseudospin representation with |ll,To) and 
| — , 5*) forming the bases. 



B. Singlet-triplet mixing 

Due to the HF interaction electrons prepared in the singlet state can be flipped to the triplet states. The spin 
state of an electron pair evolves under the Hamiltonian p6p . The initial state of nuclear spins is assumed to be in an 
ensemble, where nuclear spins are randomly oriented. Then the time evolution of the density matrix of the electron 
pair-nuclei coupled system is given as 



Pit = 0) = J2PnPn\S){S\ ^ Pit) = J2pnPn\'^nit)){^nit)\ 



|*„(t)) = (cos + sin -^)\S) - sin ^iTo) , 17„ = ^r^hl + jy2 , (38) 

where pn characterizes the nuclear spin state which assumes the HF field value: 5hzPn = hnPn, the weight of which 
is Pn, namely X^n?'" ~ Trp„ = 1. In the following h will be set to unity (/i — 1) for simplicity. From psp the 

probability to detect the triplet (singlet) state follows as 

2l2 



Pt = l( /h^ + j2 i^ Vr^h^ + JH)), (39) 

Ps = l- Prit), (40) 
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where (. . .) denotes the ensemble average over the HF fields. When the nuclear spins are unpolarized and randomly 
oriented, the spectral weight of the HF field p„ in ([55]) follows a Gaussian profile^: 

p[h] = -^e-'^'/'^' . (41) 

This is the continuum expression under the correspondence of Pn — > p[h]. — (Shi) is the mean square value of the 
HF field operator Shz- Since the nuclear spins in the left and right dots are statistically independent, we have 

(Shi) = {{hL, - hn,f) = (hlz) + {hi,) . (42) 

Thus (T^ is the sum of the mean square values of the HF fields in the left and the right QDs. 
For vanishing exchange coupling J = and r = 1, (|39|) features a Gaussian decay: 



Pt = 1/2(1 - exp[-a2iV2]) , (43) 
whereas for finite J, in the limit of i 3> J/c^, it shows a power law decaj^: 

which has been experimentally demonstrated^^. In the case of a vanishing external magnetic field, all singlet and 
triplet states are coupled via the HF interaction ((29|) . We consider the same situation, namely, an electron pair is 
initialized in the singlet state and after the HF interaction of duration t the spin state of the electron pair is measured. 
Probability for singlet detection is given as 



Ps = [{S\l/A-S,it)-S2{t)\S)). (45) 
The solution of (|29p in the Heisenberg picture yields 

Si{t) = hLhL ■ Si + (Si - fiLhL ■ Si)coshLt + Kl x Sisin/iL^ (46) 
with liL = hi /I hi I and its ensemble average over is calculated as 

(Si(0) = fl + 2(l-ai<2)exphaiiV2])Si/3. (47) 



Here = (h|^)/3 and similarly the expression for S2(t), i.e., (S2(i)) is obtained by the replacements — > and 
o'L ^ ctr- Using and (|i71) the singlet detection probability (gS]) can be readily evaluated)^ 

Ps = 1/4 + [1 + 2(1 - crit2)e-'"' _^ 2(1 - cr|ji2)e-'^«*'/2]/12 (48) 

which yields 1/3 as i ^ cx). 

When J ^ and no magnetic field is applied, the Hamiltonian is as follows: 

i/ = hi • Si + hfi • S2 + ./(Si • S2 + i) . (49) 



Within the semiclassical model, we have diagonalized (|49l) to find the singlet detection probability: 

Psit)^{\ E l(^|e,>Pe-^-*|'), (50) 



i=1...4 



where \ei), i — 1, . . . ,4: are the eigenvectors of (|49p for given and h^; values. (...) denotes ensemble averaging over 
the HF fields featuring a Gaussian distribution (|¥T|) for hi(fl) assuming aL = o-r , i.e., a = ^/ cr'i + a'j^ = V^aL(^jiy 
The time dependence of Ps{t) in ([50)1 and the asymptotic value Ps{t oo) vs J/cr are shown in Fig. [T^) and Fig. 
[T}d), respectively. The singlet probability does not feature oscillations for finite J values, in contrast to the case of 
\S) ~ \To) mixing (Eq. (HI])) which features oscillations at long time scales. This is due to the destructive interference 
between contributions from four eigenstates. 
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crt J/cr 

FIG. 1: a)Ps{t) fEa. (|50|) ) as a function of time for J/cr values of i) 0, ii) 0.5, iii) 1. h)Ps{t — > oo) as a function of J/a. These 
are calculated when no magnetic field is applied. 

C. Comparison with experimental data 

Now we make comparison of the above theoretical results with available experimental data to derive relevant 
physical parameters. In the experiments by Petta et alf^, the exchange energy J has been obtained as a function 
of the bias voltage VJ, and their experimental data are shown in Fig[2^) by dots. We fitted the experimental data to 
the expression ((27)) assuming a linear relation between the bias voltage Vf, and the detuning A/2 — qVb + Eq. The 
expression for the exchange energy becomes 

J - -Jo + qVb + Ea + VilVb + Eof + t2 (51) 

which is fitted to experimental values yielding, Jq = 5E + j — 0.07 /ieV, q/e = 5.86 x 10^^, Eq — 3.24 ^eV and 
t — 1.43 /ieV. The fitting is performed in the range A G [—9.4, l]/ieV or equivalently Vb € [— 2.15, — 0.46]mV, which is 
exhibited in Fig. [2^) with the solid line. Due to an applied magnetic field, B ~ lOOmT, the exchange energy ([STjl can 
become negative for particular values of the bias voltage. Here the singlet-triplet crossing (J = 0) occurs at a bias 
voltage Vb — — 5.53mV. 

In the experiments^^ the singlet-triplet mixing data were obtained as a function of the HF interaction period at the 
bias voltage Vb = —6 mV. For instance, the solid line in Fig. [2}3) exhibits the experimental probability to detect the 
singlet state as a function of the HF interaction period, where no external magnetic field is applied(see Eq. ([50]) ). 
In the fitting procedure first we determined the J/a ratio from the asymptotic value Ps{t ^ 1/a). In particular, 
the experimental data shown by a solid line in Fig. [^b) exhibit an asymptotic value Ps{t ^ 1/cr) ^ 0.47 which 
corresponds to the value J/a = 0.8 in Fig. [T|d). At the next step, by matching the width of the main peak of the 
Ps{t) profile for J/a — 0.8 and the experimental data, the value of cr ~ 0.09/icV is determined and thus \ J\ ~ 0.07/j,eV 
is fixed. On the other hand, we obtain J ~ — 0.038/ieV when we insert Vb — — 6mV in the expression ([5T|) . This 
discrepancy by a factor about 2 may be induced by the inaccuracy of ((5T|l because Vb = — 6mV is out of the range of 
fitting, i.e., [-2.15,-0.46]mV. However, the agreement in order of magnitude is rather satisfactory. 

III. BUNCHING OF ELECTRON SPIN MEASUREMENTS 

Now that we have formulated the basic features of the electron-nuclei coupled system, we can examine the details of 
its quantum dynamics. First of all, we reveal an interesting phenomenon of bunching of electron spin measurements 
which is caused by the correlation among successive measurements and is induced by the long-lived quantum coherence 
of nuclear spins. We also discuss the effect of relaxation of nuclear spins on this phenomenon of bunching. 

A. Successive measurements of electron spins 

Now we show that by electron spin measurements in a double QD governed by the Hamiltonian in (|36p . the coherent 
behavior of nuclear spins can be demonstrated. Corresponding to the experimentai^, we assume that an electron pair 
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-1.4 -1.2 -1 -0.8 -0.6 10 20 30 40 50 60 

H(mV) ^(ns) 

FIG. 2; a) Experimental values(dots) and theoretical fitting(solid line) for the exchange energy J as a function of the bias 
voltage Vb which is a linear function of the detuning A. b) The singlet detection probability Ps as a function of the duration 
time of the HF interaction without an external magnetic field; experiments in Ref. (solid line) and theoretical results (dashed 
line). 



is initialized in the singlet state and the nuclear spin states are initially in a mixture of Sh^ eigenstates (j38|) . In the 
unbiased regime, i.e., r = 1, the nuclear spins and the electron spins interact for a time span of r. Then the gate 
voltage is swept adiabatically, switching off the HF interaction, namely r — > 0, in a time scale much shorter than the 
duration time of the HF interaction r. Next a charge state measurement is performed which detects the singlet or 
triplet state. Probability to detect the singlet or triplet state is calculated as 

n n 

with a„ = cos il„r + iJ/2r2„ sin ri„r, /3„ = —ift,„/2il„ sin ri„T, (53) 

where the notations in l|38p are used. Subsequently one can again initialize the system in the singlet state of the 
electron pair and turn on the hyperfine interaction for a time span of r and perform the second measurement. In 
general over N times measurements, the nuclear state conditioned on k{< N) times singlet and N — k times triplet 
detection is 

'^N,k=Oj2pM'W^''-'^Pn, (54) 

n 

the trace of which yields the probabihty to have k times singlet outcomes: 

PN,k^^ra^,k^O{\anP\'^''~''^), (55) 

where (. . .) is the ensemble average over the HF field hrA- Hereafter, this case will be referred to as the coherent 
regime. One can easily contrast this result with that for the incoherent regime in which nuclear spins lose their 
coherence between the successive spin measurements and relax to the equilibrium distribution. The latter is given by 

^k.-C)(l«iy(l/3p>(^-'=). (56) 

When the nuclear spins are incoherent, the probability distribution (j56p obeys simply a Gaussian distribution with a 
mean value of fc = A^dap) and the variance of iV(|Q;p)(|/3p), as iV — > oo. However, when nuclear spins preserve their 
coherence, the probability distribution (|55p may exhibit different statistics depending on the initial nuclear state. The 
two probability distributions ([55]) and (|56p yield the same mean value, i.e., k = N{\a\'^), but with distinct higher 
order moments. If the weight factor p„ of the HF field in the equilibrium distribution has a width cr, then for the 
duration time of the HF interaction r > 1/a, the distributions ([55]) and (|56p start to deviate from each other. They 
yield the same distribution only when the initial nuclear state is in a well defined eigenstate of Sh^, i.e., when a = 0. 

If the nuclear spins are coherent over the span of the experiment, then successive electron spin measurements are 
biased to all singlet(triplet) outcomes. In particular, when the initial nuclear spins are unpolarized and randomly 
oriented, the distribution of the hyperfine field is characterized by a Gaussian distribution (|4ip with the variance a^. 
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As the simplest case, let us check the results of two measurements, each following the HF interaction of duration t. 
Probabilities in the coherent and incoherent regimes for two singlet detections are respectively calculated as 

P2,2 = = {6 + 2e-2-'*' + 8e--'*'/2}/16 

2 = (|a|2)2 ^ {4 + 8e--'*'/2 + 4e-"'*'}/16 , (57) 

where results are given particularly for J = and it turns out that P22 > P22 ■ 

In Fig. [21 for iV = 20 measurements, PN,k is shown for three values of the duration time of the HF interaction: 
err = 0.5, 1.5, 00. For t — 0, the probability for both fSS]) and ((56)) is peaked at A; = 20. However, immediately after 
the HF interaction is introduced, the probability distributions show distinct behaviors. The measurement results in 
the incoherent regime approach a Gaussian distribution. On the other hand, in the coherent case, the probabilities 
bunch at fc=0 and 20 for J = and when J/a = 0.5 those bunch at /c = 20 only. As J is increased above some critical 
value, no bunching takes place at fc = times singlet measurement, since the singlet state becomes energetically stable 
and the state change to the triplet state becomes unfavorable. 

To observe the bunching N successive spin measurements are performed within the coherence time of the nuclear 
spins. Then after waiting for some time so that nuclear spins are again randomized, another set of N successive 
measurements are carried out and so on. Thus an ensemble average of N measurements is performed which results in 
a bunching of either spin singlet or triplet outcomes. This bunching is a clear signature of coherent behavior of nuclear 
spins, which can easily be contrasted with the incoherent regime which merely exhibits a Gaussian distribution. 



B. Effects of nuclear spin diffusion 



Now we will discuss the effect of nuclear spin diffusion on the bunching of electron spin measurements, which 
leads to a transition from the coherent regime to the incoherent regime. During the interval between the successive 
measurements the nuclear spin state relaxes to the equilibrium distribution due to the dipole-dipole interactions^^. 

If the substrate surrounding a QD is of the same kind of material as that of the QD, nuclear states will diffuse due to 
the interaction with the surrounding nuclei, which leads to a change both in the total spin angular momentum of the 
nuclei and the HF field. The inhomogeneous distribution of HF coupling constants will also induce a redistribution 
of the spin angular momentum leading to a change in the HF field. Since a detailed discussion on the nuclear spin 
diffusion is beyond the scope of this paper, we simply develop a phenomenological argument based on the diffusion 
equation in the phase space of the HF field: 

^^ = ^^ip[h,t]-po[h]). (58) 
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FIG. 4: Results of two successive spin measurements: Probability of a) two triplet b)one singlet-one triplet measurements, at 
=i)0, ii)0.01, iii)0.05, iv)0.1, v)0.5, and vi) oo. Inset: (Measurement scheme) r is the duration of the HF interactions, 
and T the waiting time with HF interaction switched off, where electron spin measurements are denoted by M. 



where po[h] is the distribution corresponding to the steady state configuration of nuclear spins. At high temperatures 
compared with the nuclear Zeeman splitting, po[h] obeys a Gaussian distribution (Hi}. The general solution of the 
diff'usion equation ([55)) can be cast into the form: 



p[h, t] = _L=e-i^ - , ^ e ^i^^ + -yL= I dh'e-'-^p[h', t = 0] , (59) 

V27rcr2 ^27r(CT2 + 2Kt) ^A-kkI J 

where p[h,t = 0] is the initial distribution of the HF field. 

The randomization of nuclear spins will lead to loss of memory effects described in the last section. The nuclear 
state conditioned on the electron spin measurements (Eq. (I54p ) will decohere throughout the successive measurements. 
The nuclear diffusion time is much longer than the characteristic time of the spin singlet-triplet mixing induced by 
the HF interaction (Eq. ([36]) ). namely tdif f t ^ l/5hz- Decoherence of nuclear spins will mainly take place during 
the electron spin measurement, because this process of spin-charge conversion is time-consumingl2. 

For instance, when an electron pair is initialized in the spin singlet state starting with a randomized nuclear spin 
configuration (PT|) . then subject to the hypcrfinc interaction ([55)1 for J = 0, of duration r and is followed by the 
electron spin measurement, the spectrum of the HF field becomes 

p[h] = A/'e-'*'/2<T^ (1 ± cos hr) , (60) 

corresponding to either singlct(+) or triplet(-) outcome, where A/" is a normalization constant. Governed by the 
diffusion equation ([55| , the distribution (pO)) after a time span of T evolves to 



1 ^2 p 2(2»T + ,t2) ^ 1 ± e 2»T + <t2 cos ..J^Tp, -i 

p[h,T] = ^=e-t^ + ; — I l+2.T,a^ 

VW ^27r(2AtT + a2) I l±e^ ^ 

This distribution converges to a Gaussian for T 3> a'^ / n. It is also to be noted that the duration of the HF interaction 
T also affects the effective diffusion time. This means that the period of modulation ~ 1/t induced in the nuclear 
field spectrum affects the speed of diffusion. In fact (|6T|) approaches the Gaussian form pQ[h] when t'^ ^ l/'«7', l/f^, 
namely when the period of undulation in the nuclear field spectrum is short enough to be smoothed out easily, the 
distribution converges to po [h] ■ 

After a time span of T following the first measurement, the system is again initialized, HF interaction is switched on 
for a time span t', then a second spin measurement is performed. Here typically T 3> t, r'. The measurement results 
approach that of semiclassical picture when T ^ ct^/k or ^ ^/kT, Otherwise one can still trace the nuclear 

memory effects in the measurement results. In Fig. |4l some examples are shown for two successive measurements with 
parameter values of kT/ct^ — 0, 0.05, 0.1, 0.5, oo, as a function of the HF interaction time r = t'. In the asymptotic 
limit of kT/ct^ = oo, we can check that the probability in Fig. [l]a)(b)) approaches -P2 o(A' i)- 
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FIG. 5: Measurement scheme to observe the electron spin revivah Each time the electron spin is initialized in the spin singlet 
state (denoted by S), then the HF interaction is switched on for a period Ti followed by the electron spin measurement with an 
outcome Mi, for i = 1, • ■ • , A'^. These are the preparation stage. Then after the HF interaction for a period t, the + 1-th 
measurement is carried out. 



IV. PURIFICATION OF NUCLEAR SPIN STATE AND ELECTRON SPIN REVIVALS 

In this section we investigate the conditional preparation and purification of nuclear spin state via successive electron 
spin measurements. This feature becomes manifest via revival phenomena of the electron spin state. Here the HF 
interaction is assumed to take place in the unbiased regime of the double QD, i.e., when J — and r = 1 in (|38p . 
Then the nuclear state prepared by N successive electron spin measurements with k times singlet outcomes, each 
following the HF interaction of duration times ti,T2, . . . ,tn is given by 



hnT, 



N 

n ^i'^ 

j=k+i 



2 ^n^j 



(62) 



where A/" is a normalization constant. The sequence of measurements is depicted in Fig. [5l In the following we 
consider the case where all measurement outcomes are singlets and examine two typical cases: A. ti = r2 = . . . = tjv 
and B. n = 2t2 = . . . = 



First case: ti = T2 — ■ . ■ — tn ~ t 



In this case all the duration times of the HF interaction are equal and the prepared nuclear state following N times 
singlet measurements is given by 



a PnPn COS — - 



(63) 



where A/" is a normalization constant. Given the initial state p{t = 0) = <t|S')(S'|, the probability to measure the 
singlet electron spin state after the HF interaction of duration time t is calculated as 



P{t;{n = t}^=i,..., 



N 



(cos^^ [fer/2] cos'[ht/2]) 
(cos2^[/ir/2]) 

1 l E'foE.^±ff)(exp[^(.~jV)fer + ^afet]) 

2 4 E^fo(?)(expW^-A^)H> 



(64) 



where (. . .) denotes ensemble averaging with respect to the initially random nuclear spin state (|4ip . Using the identity: 

(e^"*) ^ [ dfte-'^'/^^'e*"* = e-'^'*'/^ , (65) 

the equation ([M)) can be cast into the form: 



1 , lE^" 



2N (2N\ -a^{t-{N-s)Tf/2 



P(i;{T, = T},=i,...,jv) = 7r + :r 



s=0 \ s 



2 2 ^^^Q (■2Af-Jg-(s-Ar)2<T2T2/2 



(66) 



For i < 1/(7 this gives a Gaussian decay(see Fig. |6^), whereas for t > l/a it exhibits revivals(see Fig. [B}:). For ar ^ 1 
the expression ((66l) reduces to 



^l,...,N) 



1 1 Es^o (2^)e-'"'(*-(^-")^)'/2 

2 ^ 2 



(67) 
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FIG. 6: Conditional probability for singlet state detection as a function of HF interaction period at, subject to A'' = 0, 1, 2, 5, 10 
times prior singlet state measurements and for HF interaction duration times a) ar = 1.0, b) err = 3.0, and c) ar = 6.0. 



featuring revivals at t = nT{n = 1,2,...) with a decreasing amplitude: 

'/-(;f„)K«) 

which becomes 1/2 + e~" /^/2 for N ^ 1. In the method proposed here, the nuclear spin state can be conditionally 
purified without determining the precise value of the HF field. Although the HF field may be still assuming indefinite 
values, electron-nuclei correlations lead to revivals at known times. As an example consider the case when the nuclear 
spin state is prepared by five HF interaction stages, each of which has a duration time r = lO/cr and is followed by 
a singlet detection of the electron spin state. This conditionally prepared nuclear spin state revives the spin singlet 
electron state at times at — 10, 20, . . . , 50 with fidelities 1/2 + (^'^^) /2(^j^) , for = 5 and s ~ 1, 2, . . . , 5 which are 
11/12, 31/42, . . . , 253/504. Success probability to prepare such a state is ~ 1/2^. 

In order to understand the physics of the revival more clearly, we consider the limit A ^ 1. Since | cos 6*1 < 1, 
cos^^ 9 is sharply peaked at 6 — sir (s G Z) and can be approximated as 



cos^^ 9 



= ^(1 - N{B - sirf + ---)c^Y^ exp[-A(e' - sT:f] . (69) 

seZ sGZ 

Then the spectrum of the nuclear HF field corresponding to the conditionally prepared state (|63p can be approximated 
in the limit A^ ^ 1 as 

p[h] = AAe-'''/2-' cos^^ — ~ AAe-'''/2.^ ^ e-^'^-h.^V^.l (70) 

with = , (71) 

implying squeezing of the HF spectrum at particular known values hs = 2stt/t. Given the initial nuclear spin state 
with the spectrum ((70|) . the probability to recover an initial singlet electron spin state after the HF interaction of 
duration time t is given by 

P{t) = (cos^ ht/2) = l + ^J2 e-'^"*'/2e""'/2-^ cos hst , (72) 

where it is assumed that a ^ am = l/(T-y/A^/2) and the normalization constant Af is set to satisfy P{t = 0) = 1. 
In (I72p each gives rise to revivals at times t^ = t/2s and its integer multiples. At the common multiples of all t^ 
values which are nr, n — 1,2,3..., the probabilities add up coherently leading to revivals(c.f. ()66|) ) and each revival 
has the amplitude 1/2 + e-"'/^/2. 

The revival phenomenon also applies to some arbitrary initial electron spin state subject to the HF interaction with 
the conditionally prepared nuclear spin state ((63|l . When the initial state of the system is assumed as 

p(i = 0)=a|V)(^|, 1^) =cos^|^)+sin^e-*^|T) , (73) 
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the fidelity F = {ilj\p{t)\ilj) to recover the initial electron spin state \ip) at time t is calculated as 

F = sin^ 9 cos^ + (1 - sin^ 9 cos^ (t>)Pit) (74) 

with P{t) given by (|64l) . 



B. Second case: ri = 2r2 = . . . = 2^ ^rjv = r 

Here the duration times of the HF interaction are decreased by one half successively, namely t = ti — 2t2 = 2^x3 = 
. . . — 2^^^tn- In this case the prepared state (p^ is given as 

a=MY,PnPn'[[cos'^ (75) 

n s— 1 

and the probability to recover the initial singlet electron spin state after the HF interaction of duration time t is 
calculated as 

The singlet state is revived a.t t — 2t x Q.I1I2 ■ ■ - In, where k ~ 0,1. This amounts to 2^ revivals at times t = 
t/2^-1, 2r/2^-i, 3r/2^-i, . . . , 2r(l - 1/2^), 2r. 

Now we briefly discuss the HF spectrum of the state ([75]) and its relation to revivals in ([75]) . The HF spectrum of 
([75)1 can be cast into the form: 

^ bill 2N 

In the same way as in (iTOl) . in the hmit of A'' 1, we can show that 

p[h] = AAe-"'/2'^' e-("-''=)'/'"" (78) 
sez 

.AT., , /3 1 



with hg = s2 tt/t and dm = A/ — • (79) 

V 2 T 

This implies squeezing of the spectrum at hg — s2^tt/t, s G Z. Given the initial nuclear spin state with the spectrum 
([78)) . the probability to recover the initial singlet electron spin state at time t is given by the same expression 
as in ([7^ . Each hg leads to revivals at times r/(s2^) and their integer multiples, which add up coherently at 
nT/2^^^, {n — 1,2, ■ ■ ■), giving rise to revivals (c.f. (j76p ). Thus we can understand that the revival phenomena occur 
reflecting the undulation in the nuclear field spectrum induced by the electron spin measurements. 

As a concrete example we make a comparison of the two schemes for N = 2, k = 2 and examine the electron spin 
revivals for this conditionally prepared state. For ti = 2t2 — t ^ l/f, the conditional probability (|76p is given as 

P{t;{Ti=T,T2 = -})-- + -{e -2 +2e ^ + 3e 2 + 4e -} , (80) 

whereas for n = r2 = r ^ l/cr, (|66l) is calculated as 

P{t; {n = r, r2 = r}) ~ - + -{e ^ + 4e — + 6e-— } . (81) 

We have more revivals with higher probabilities for the former case in which the undulation in the nuclear field 
spectrum is more structured. 

In the above we found that the nuclear field spectrum is squeezed or undulated through the electron spin mea- 
surements. In order to examine the degree of squeezing quantitatively, we estimate the purity of the nuclear spin 
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FIG. 7; (Color online) Purity (in units of I/O) as a function of the duration times ri and T2 of the HF interaction in the case 
of A'^ = 2 measurements for a)fe = 2 and h)k = times singlet outcomes. 



system. The purity of the system is given by 'PN,k = Tr cr^ f.. Using the identity Tr pf^ = l/Vpn, where V is the total 
dimension of the Hilbert space of nuclear spins, we obtain the purity of the state ([62|l as 



'''' ^(/d/.M^]ntrSin^^nf..+i-n^^)^ 



Sfc+1 + ...+sjv 



1 St^oU)U)---(/Je-^['-^~^)^^+(--^)^-+-+<--^)^"]^~l) 



(83) 



where f = err and sums were evaluated in the continuum limit X)nP" ^ /d/ip[/i] in (|82p . We can extremize the 
Duritv(|83p by choosing appropriate duration times of the HF interaction . In the asymtotic limit 1, we have 

^ 1 EUo O ■ ■ ■ - 2)n + . . . + js^ 2)r^] 

' ^ E'.=o O ■ ■ ■ (/J(-l)^^-^+-+^"'5[(si - i)n + . . . + (.A. - 1)%]]' 

From (|84p we see that there are several asymptotic values determined by the roots of the linear equations ^(s — 2)Ti = 
and J2is^ - = 0. 

For instance in the case of = 2 and A; = 2, we sec that the numerator in ([M]) has a contribution from the choice 
of si — S2 = 2, whereas the denominator has a contribution from the choice of si = S2 — 1, irrespective of the relative 
magnitude of ri and T2. Additionally, in the case of ri = 2t2 or T2 = 2ti the numerator has finite contributions 
arising from some combinations of si and §2. On the other hand, the denominator does not have such contributions, 
because the equation '^^{si — l)Ti = cannot be satisfied except for si = S2 = 1. In the case of ri = T2, both the 
numerator and the denominator have finite contributions from appropriate choices of si and S2 other than the trivial 
ones given by si = S2 = 2 or si = S2 = 1. Summarizing, there are three asymptotic limits (see Fig. [7^ ), namely 
when i) ri — 2t2 then 7^2,2 = 11/41?, ii) ri — T2 then 7^2,2 — 35/181? and c)otherwise 7^2.2 = 9/41?. In general the 
purity attains its maximum for all singlet outcomes, i.e., for k ^ N and under the condition that the duration times 
of the HF interaction are halved at each step, viz., Ti/ri^i — 2. 

In Fig. [5] the purity Vn.n is shown as a function of the number of measurements N , in the asymptotic limit of 
Ti ^ 1/cr, i = 1, . . . , iV. The curve i) corresponds to the maximum purity and the curve iv) to the minimum, whereas 
all other choices of interaction periods ti : T2 tj^ yield intermediate values (see Appendix [C]). 
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FIG. 8: Purity Pjv,jv(in units of 1/T>) is shown in the natural logarithmic scale, for the schemes i)ri = 2t2 — . . . — 2^~^tn, 
ii)ri = 2t2 = 3r3 = . . . = Ntn, iii)Ti : tj is irrational for any pair of and for iv) n = r2 = . . . = rjv as a function of the 

number of electron spin measurements(A'' = 1, 2, . . . , 8). 

V. REALIZATION ON A SINGLE QD 
A. Single electron on a single QD 

So far we have discussed the bunching and revival phenomena only for a double QD system. The same predictions 
can be made also for a single QD occupied by a single electron22i21i2^. Consider a single QD occupied by a single 
electron, under an external magnetic field such that the electron Zeeman energy is much greater than the HF energies. 
Then the system is described by the Hamiltonian: 

H ~ gefiBBS, + h,S, = [B, + K)S, , (85) 

where gg, is the electron g-factor, the Bohr magneton, B the external magnetic field applied in the z direction, B,, 
represents the electron Zeeman energy and is the nuclear HF field in the z direction . Spin fiips are suppressed 
since B^ = gefJ-sB ^ -y/ (h^). The spin eigenstates in the x direction |±) = (| t) ± | J,))/\/2 are coupled by the HF 
interaction with | f (J,)) being the eigenstates of 5*2. The time evolution of the state |+) is simply given by {h — 1) 

e-^^* 1+) = cos ^^±^t\+) - i sin ^^±I^t\-) . (86) 

Now let us consider the following experiment. Each time the electron is prepared in the state |+). Next it is loaded 
into the QD, then removed from the QD after some dwelling time r and the spin measurement is performed in the 
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basis of |±). Essentially the same predictions as those for a double QD can be made for this system, namely the 
electron spin bunching and revival. We consider the electron spin revival as an example. After N times the HF 
interaction of duration time r, each followed by the measurement outcome of the |+) state, the nuclear spin state 
becomes 

E2N V^e + hnjT . 
PnPn coa^ , (87) 



where the initial distribution p„ characterizes the random distribution of the HF field (|4T|) . Then the electron spin 
state is prepared in 1+), yielding the initial state p{t = 0) = which evolves under the Hamiltonian ([55]) . The 

probability for obtaining 1+) after the HF interaction of duration t is given as 



(C0S2^(^4^C0S2 (^4^) 



(C0S2^ (M)r ^ 



1 , iE..±E?fo(?)(^^^"=^"^'^ 



s-N)T+at]\ 



2 4 ^2Af^((2Ar^gj(s_Ar)(B,+/i)r^ 



where (. . .) denotes ensemble averaging with respect to (|4T|) . Using the identity ((65|) the equation ([88|) for t ^ 1/a 
can be put in the form: 



2N 



Pit; {n = rh=i,...,N) 0, 1/2 + 1 f](2f )e-^(*-(^-)-)V2 cos B,[t - (N - s)r] . (89) 
Thus, in a single QD, revivals are present as in the double QD case(c.f. Eqs. and (|67p). 



B. A pair of electrons on a single QD 

The Hamiltonian ((36)) can also be used to describe a pair of electrons in a single QD^Si^S and the same predictions 
as those for a double QD can be made. In the two electron regime, the energy splitting between the singlet ground 
state and the triplet excited state can be tuned down to zero by application of a magnetic field^i^i^ leading to a 
singlet-triplet crossing. Under a high magnetic field, the triplet state (Tq) having zero magnetic quantum number is 
coupled to the singlet state (S) via the HF field: 

h = Av„Y,M^^)€i^^)Ii'^/^: (90) 

i 

where (t>g[e) is the ground (excited) state orbital in the QD and the derivation is given in Appendix A. Typically for a 
two-dimensional QD with harmonic confinement, the HF field (j90p has a mean square value: 

{h^fh) = A^I{I + l)vo/lQndrl , (91) 

where tq = ^Jh/mVL is the Fock-Darwin radius, d the thickness of the QD, VL = ^Juj^+uj^JA with LdQ{ujc) being the 
frequency of the harmonic confinement potential(the cyclotron frequency) and / is the magnitude of the nuclear spin. 
In the energy spectrum of a single QD the singlet-triplet crossing was observed via the tuning of magnetic field^. 

For an isotropic GaAs QD with a harmonic confinement energy uj = ImeV, a singlet-triplet crossing will take place 
at S ~ I.IT and the second excited state which is a single t, is se parated by ^ 0.2meV. For such a QD, with thickness 
d = Inm, the rms value for the HF field ([91]) will be y/JJiXh) ~ 0.04 /icV, which implies that the system can be 
treated as a two level system coupled by the HF field(See Appendix |B] for the spectrum of a single QD occupied 
by two electrons). The relevant Hamiltonian describing the dynamics within the subspace formed by \S) and |To) 
is essentially the same as that for an electron pair in a double QD. Furthermore, the electrons' spin state can be 
initialized and measured with high fidelity by a spin-selective coupling to leads, relying on the spin-dependent tunnel 
rates^. Thus the observation of the same phenomena as the bunching in the electron spin measurements and the 
revival of the initial electron state is feasible also in a single QD occupied by a pair of electrons. 
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VI. CONCLUSION 



We have investigated the quantum dynamics of the electron-nuclei coupled spin system in QDs and predicted some 
interesting new phenomena. The quantum correlation induced in the system via consecutive HF interactions leads to 
the bunching of outcomes in the electron spin measurements and the revival of an arbitrary initial electron spin state. 
Simultaneously, the nuclear spin system is affected by the quantum correlation and is in fact squeezed as confirmed 
by the increase in the purity. It is suggested that the consecutive electron spin measurements provide a probabilistic 
method to squeeze or prepare the nuclear spin system. We also discussed the effect of nuclear spin relaxation on the 
bunching and revival phenomena based on a phenomenological model and exemplified a change from the coherent 
regime to the incoherent regime. All the results obtained are applicable not only to a double QDs occupied by a pair 
of electrons but also to a single QD occupied by a single electron or a pair of electrons, whenever the HF interaction 
is present and the nuclear spin state is coherent throughout the experiments. 
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APPENDIX A: HYPERFINE INTERACTION FOR AN ELECTRON PAIR IN A SINGLE QD 

Here we derive the Hamiltonian for an electron pair in a single QD. Under a sufficiently strong magnetic field, 
the triplet states T± are well separated from the Tq state and the singlet state S. Thus the Hamiltonian within 
the subspace spanned by Tq and S states will be considered. The wavefunctions for the S and Tq states are given, 
respectively, as 

*s(ri,Ci,r2,C2)=^,(ri)03(r2)-i=(a(a)/3(6)-/3(a)«(6)), (Al) 
vI/To(ri,a,r2,6) = ^{M^i)Mr2) - Mri)M^2)){a{ii)m2) + miM^2)) , (A2) 

where (j)g{(j)e) is the ground(excited) state orbital in the QD and a(/3) denotes the spin up(down) state. The HF 
interaction for two electrons is given by 

Vhf = Avo J2 Si ■ I< -^(ri - + Avo S2 • I. 5(r2 - R,) , (A3) 

i i 

where Ri(I^) denotes the position(spin vector) of a nucleus and Si and S2 are the electron spin vectors. Then we find 

{^To\Vhf\^s) = -^AvoJ2€(r^)M^'^)I^z , {^s\Vhf\^To) = {^ToIVhfI^s)* , (A4) 
^ i 

i^slVHFl'i's) - (*t„|Vhf|*To> = . (A5) 

Thus the singlet-triplet mixing is induced by the HF interaction. The effective nuclear field operator coupling the 
singlet and triplet states in (jASp will be introduced by 



h^{^To\VHF\'^s) (A6) 

which has the dimension of energy and its mean square value is estimated as 

(/.ftt) = i^j2mR.)^A^^nii) (A8) 

i 

A^voI{I + l) /rf3^|^*(r)^^(r)P, (A9) 
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where / is the magnitude of the nuclear spin. Employing the envelope functions for the ground and excited states 
given by 

(l)g{r,e,z)^-j=-e ^^-cos( — ), (AlO) 
Ur.e.z) = ^e"4 re-*^/Icos(^) (All) 



with r. = xl^,^-p.^i^?. (A12) 



where d is the thickness of the QD, we have 



m-^.m^l). (A13) 



APPENDIX B: ENERGY SPECTRUM OF TWO ELECTRONS IN A QD WITH ISOTROPIC 

HARMONIC CONFINEMENT 

Here we calculate the energy spectrum of two electrons in a QD, assuming an isotropic harmonic confinement of 
frequency in the xy plane and a strong confinement along the growth(z) direction. Introducing the center-of-mass 
and the relative coordinates, the total Hamiltonian can be divided as 

H = HcM + H,ci. + Hz , (Bl) 

Hz=g^iBB{SP +SP) (B4) 
with R= (ri +r2)/2, P = pi +P2 , r = ri - ra , p = (pi - P2)/2, (B5) 



^(CM) ^ R X P , L('-<=') = r X p , n= ^ujI+ujI/A , (B6) 

where ri(2) denotes the coordinates of the first(second) electron in the xy plane, M — 2m*, /i — m*/2 with m* 
being the electron effective mass, lOc = eB/m*c the cyclotron frequency and is the Bohr magneton. Employing 
TO* = 0.067 • TOe, K = 12.53 and g — —0.44^^ appropriate for GaAs, we diagonalized numerically the Hamiltonian 
Hici. for the relative coordinate part. In the numerical diagonalization, 20 Fock-Darwin basis functions are employed 
to guarantee sufficient accuracy. The energy spectrum for a GaAs-like QD with ujq = ImeV is depicted in Fig. [5^) 
as a function of the magnetic field B for the orbital part Hqm + -ffrei- In Fig. the energy spectrum is plotted 
in the vicinity of the lowest energy singlet-triplet crossing point including the spin degrees of freedom. The singlet 
ground state and the triplet first excited states feature a crossing at _B ~ IT separated from the next excited state 
by ^ 0.2meV which is a singlet. For a magnetic field B < 0.9T, the electrons can be loaded into the singlet ground 
state and then by sweeping the magnetic field to the S — Tq crossing point, the system can be initialized. Here the 
S — r_(- crossing point should be passed at a rate much faster than the HE interaction time which was estimated as 
~ 103 ns from ^/JJihX) ^ 0.04/ieV in Sec. V-B. At the S* — To crossing point the T_|_ and r_ states are separated by 
a Zeeman energy ^ 25^eV which is much greater than the HE interaction energy. Thus at the S — Tq crossing point 
the system can be described by a two level Hamiltonian composed of \S) and jTo) states. As demonstrated in the 
experiments by Meunier et al.— , the phonon mediated spin relaxation time exceeds well beyond ms order which leaves 
the HE interaction as the only relevant mechanism at time scales shorter than ms. 



APPENDIX C: PURITY OF THE NUCLEAR SPIN STATE AS A FUNCTION OF THE NUMBER OF 

ELECTRON SPIN MEASUREMENTS 

When all : Tj, (i 7^ j) ratios are irrational, in (|84p. only si = S2 = . . . = sjv = 2 in the numerator and only 
si = S2 = . . . = Sat = 1 in the denominator contribute, yielding Vn.n — (3/2)^/1? (Eig[Hl iii). 
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FIG. 9: a)Energy spectrum of two electrons in a GaAs-like QD with a harmonic confinement energy ujo = ImeV. Sohd(dashed) 
hnes indicate the triplet (singlet) states, b) Energy spectrum in the vicinity of the singlet-triplet crossing. Triplet states T± are 
split by the Zeeman energy ±gHB — =F25.5/^eV/T from the To state as shown by solid lines. 



In case of ri = T2 = ... = tat = r. the purity is given by Vn,n ~ {2%) I Cn) /-^ (FiglHl iv). This can be verified by 
inserting Vr^ = r in 7'jv.jv(Eq. (|82p) then taking the limit t ^ 00, 
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where in the last line hmT-^oo(exp (ihriT)) = hnv_>oo exp (— n^cr^T^/2) = Sn.o from (j65[) is employed. 

For the case ti = 2t2 = 2^X3 = . . . = 2^~^T]\[ = t, the asymptotic value of the purity in the limit of err —> 00 can 
be evaluated more systematically from the expression in (76) rather than (82). The density matrix after the N times 
measurements is given by 
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where the normalization constant M is determined by 



1 = Tr pN,N = TN^Pn- 
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The last factor takes a large value about 4 near /i„ ^ hs — 2 stt/t{s G Z) and can be approximated as 
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In the limit of (tt ^ 00, the Gaussian distribution pn is much broader than the last factor in (jC4p and we have 
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Assuming furthermore (jt ^ 2^ , the sum over the integer s can be replaced by an integral and Af can be fixed as 
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This result is equal to the exact result Af = 2 , i.e., 
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Now that the density matrix is determined, the purity is calculated as 



Vn.n = Tr CT^^jv ^^T^J2p^ 



Vm ^^"sin4(/i„r/2^) 



By the same arguments as above, we can approximate the last factor as 
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and under the condition ut ^ 1 we have 
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Assuming err 3> 2 , the summation over s is replaced by an integral and we obtain 
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V 16^ ^''"sin4(/i„r/2^) I? 3 
This expression reproduces very well the result in FiglHl i. 
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